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SUMMARY 


bloud  physics  measurements  by  the  Twin  Otter  aircraft  are  affected 
to  some  degree  by  local  flow  effects  from  the  wing  and  from  the  probe 
canister.  Droplet  images  and  concentration  are  distorted  in  the  worst  case 
by  no  more  than  20%  for  an  aircraft  lift  coefficient  (CL)  of  0.75  and  by  only 
a  few  percent  for  low  values  of  CL.  Corrections  for  a  Particle  Measuring 
Systems  probe  mounted  on  a  short  pylon  at  the  inboard  span  wise  location 
have  been  given  which  are  based  on  an  experimentally  verified  model  for  the 
flow  velocity. 


SOMMAIRE 


Quelques  donnies  physiques  des  nuages  mesurees  par  l’avion 
‘Twin  Otter”  sont  modifies,  a  un  certain  degr6,  par  les  effets  d’6coulement 
d’air  local  autour  de  l’aile  et  de  la  boite  contenant  la  sonde.  Les  images  et  les 
concentrations  mesurees  des  gouttelettes  sont  diformfes  au  phis  de  20% 
pour  un  coefficient  de  portance  (CL )  de  0.75  et  de  quelques  pour-cent  pour 
les  petites  valeurs  du  CL. 

Les  corrections  qui  sont  pr£aent£es  pour  les  donnfes  ont  et4  obte- 
nues  avec  la  sonde  “Particle  Measuring  System”  mont£e  sur  un  pylon  court 
fix5  a  un  point  interne  de  l’aile,  ces  corrections  fluent  dfrivges  d’un  modele 
pour  la  vitease  de  l’4coulement  et  v&rifttes  exp^rimentalement. 
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AIRCRAFT  FLOW  EFFECTS  ON  CLOUD  DROPLET  IMAGES  AND  CONCENTRATIONS 


1.0  INTRODUCTION 

The  Twin  Otter  aircraft  has  been  instrumented  as  a  cloud  physics  research  aircraft  (Ref.  1) 
and  has  been  used  for  a  number  of  years  in  co-operation  with  the  Atmospheric  Environment  Service. 
The  aircraft  can  measure  atmospheric  state  variables  and  air  motion  in  real  time.  In  addition,  cloud 
droplets  and  particles  can  be  sized  and  imaged  by  Particle  Measuring  Systems  (PMS)  probes  that  are 
mounted  under  the  wings.  It  was  known  that  the  measurements  by  the  probes  could  be  influenced  by 
the  flow  around  the  wings  and  by  the  cylindrical  body  (canister)  which  houses  the  working  compo¬ 
nents  of  the  probe.  The  aircraft  and  PMS  probes  are  shown  in  Figures  1  and  2. 

A  theoretical  model  of  local  flow  velocity  at  the  sampling  point  of  a  PMS  probe  has  been 
developed  (Ref.  2)  and  experimentally  verified  by  flight  measurements  using  a  pitot-static  equipped 
PMS  canister  mounted  under  the  Twin  Otter  wing  at  four  different  positions.  This  model  will  be  used 
in  assessing  local  flow  effects.  It  is  felt  that  meaningful  conclusions  on  cloud  physics  measurements 
can  be  made  as  a  result  of  this  experimental  verification. 

This  report  discusses  implications  of  local  flow  effects  on  cloud  physics  measurements 
related  to  droplets.  Results  for  non-spherical  particles  cannot  be  accurately  calculated  because  such 
particles  can  be  widely  variable  in  shape  and  can  have  an  unpredictable  angular  orientation  while 
passing  through  the  probe  sample  volume.  Further,  the  density  of  ice  particles  or  snowflakes  can  vary 
from  0.2  gms/cc  (say)  to  about  1  gm/cc.  The  motion  of  a  particle  is  determined  by  its  drag-to-mass 
ratio;  hence,  variations  in  mass  and  aerodynamic  properties  of  non-spherical  particles  make  accurate 
calculations  of  their  motion  impossible. 

In  the  sections  that  follow,  an  analysis  of  aircraft-induced  errors  in  cloud  physics  measure¬ 
ments  will  be  made  for  water  drops.  We  will  first  develop  the  equations  of  motion  for  the  drop,  pose 
the  boundary  value  problem  for  drops  being  sampled,  solve  the  problem  numerically  and  present  re¬ 
sults  for  drop  velocity  and  trajectory  expansion  at  one  probe  location.  A  method  to  correct  cloud 
physics  results  obtained  from  the  Twin  Otter  will  be  presented. 


2.0  EQUATIONS  OF  MOTION 

The  co-ordinate  system  (Fig.  3)  is  fixed  to  the  mid-chord  point  of  the  aircraft  wing,  x'  and 
i  being  orthogonal.  The  x'  axis  is  the  line  joining  the  leading  and  trailing  edges  of  the  wing  and  it  is 
positive  towards  the  trailing  edge.  This  axis  is  inclined  at  an  angle  a  to  the  flow  at  infinity.  The  z  axis 
is  positive  up.  The  droplet  has  velocity  components  U,  W  in  the  x',  z  directions  and  the  local  flow 
velocity  has  components  U,,  Wa  with  respect  to  these  axes.  The  equations  of  motion  for  the  drop  in 
the  x',  z'  axes  are 


m 


dU 

dt 


-  D, 


dW 

m—  -  -  Dz<  -  mg 


(1) 


The  motion  of  the  drop  is  considered  to  be  two-dimensional,  the  span  wise  positional 
variation  being  neglected  in  the  present  case.  For  probe  locations  near  the  wing  tip,  this  restriction 
would  not  be  applicable.  The  flow  about  the  Twin  Otter  wing  is  virtually  two-dimensional  at  the 
present  PMS  probe  span  wise  mounting  locations  (0.618  and  0.747  of  semispan). 


The  x'  and  i  components  of  the  drag  are  written  in  terms  of  the  droplet  velocity  with 
respect  to  the  air  in  the  form: 


V-.1 


Dx.  -  D 


(u-u,) 

/( u-  UJ2  +  (W-  wa)2 


Dz.  -  D 


(W- wt) 

/(U-Ut)2  +  (W- Wa)2 


Following  Reed  (Ref.  3),  the  drag  (D)  is  written  in  terms  of  the  drag  in  Stokes  flow  (Df). 


D  =  D(D,/Dt) 


-  (CD/CDj)3irM5  V(U-  U,)r  +  (W-  Wg)2 


where  CD  /CD  is  a  function  of  the  Reynolds  Number  and  8  is  the  drop  diameter.  The  absolute  air 
viscosity  is  n.  The  Reynolds  Number  is 


Re  -  (5/v)  V( U-  U,)2  +  (W-  W,)2 


where  v  is  the  air  kinematic  viscosity. 

Substituting  (3)  and  (2)  in  (1)  and  writing  the  drop  mass  in  terms  of  drop  diameter  and 
density  (pd )  leads  to: 


dU  /  CD  \  18  p 
dt  \CD./  52  Pi  ^ "U* 


dW  /  CD  \  18  p 

—  —  —  KW-  W  ) 

dt  \CD./  82  Pi 
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Equation*  (5)  and  (6)  are  best  handled  in  a  non-dimensional  form  where  velocities  are 
measured  relative  to  the  air  velocity  at  infinity  (V),  lengths  are  measured  relative  to  the  quarter-chord 
length  of  the  wing  (a)  and  a  non-dimensional  time  r  is  defined: 


u  »  U/V 

w  -  W/V 

x  =  x'/a 

z  ■  z'/a 

r  -  t/t* 

t*  -  a/V 


Substitution  of  (7)  in  (5)  and  (6)  leads  to 


dx 

dr 


*  u 


(7) 


(8) 


da 

dr 


(9) 


du 

dr 


~  C(Re)  K,  (u-u.) 


(10) 


dw 

dr 


-  C(Re)  K,  (w-  wp)  -  K2 


(ID 


where  the  following  non-dimensional  quantities  have  been  defined: 


C(Rc)  =  CD/CI)t 


(12) 


18  p  p 

(5/a)2  Pa  Va 


(13) 


K2  -  ga/V2 


(14) 


Re  *  ~  (5/a)>/(«'  ua)2  +  (w-  wg)2 


(15) 


u(0)  «  cos  a 


w(0)  =  sin  a 


The  trajectory  from  initial  to  final  position  is  obtained  by  numerically  integrating  Equations 
(8)  to  (11)  using  the  data  from  Reed  and  Vennard  (Table  1)  to  interpolate  for  C  and  a  Runge-Kutta- 
Gill  integration  routine  with  variable  step  size.  The  standard  atmosphere  is  assumed  (Ref.  4,  p.  21) 
where  temperature,  density  and  static  pressure  are  determined  by  altitude  h  relative  to  the  sea  level 
(SL)  values: 


Tsl  -  3.566  *  (h/1000) 


p  =  pSL  exp(-  0.0297  *  h/1000) 


exp(-  0.0347  *  h/1000) 


The  variation  of  kinematic  viscosity  with  altitude  is  obtained  to  a  very  close  approximation 
by  linear  interpolation  between  the  entries  in  Table  2.  This  table  also  contains  the  sea  level  values  of 
the  variables  in  Equations  (17)  to  (20). 

The  boundary  value  problem  is  solved  by  a  "shooting  method”  where  the  initial  position  of 
the  drop  is  chosen  and  the  final  position  is  calculated.  If  this  final  position  is  not  sufficiently  dose  to 
the  desired  final  position,  then  the  initial  position  is  chosen  again  and  the  procedure  is  repeated.  A 
simple  algorithm  for  choosing  the  initial  position  was  developed  by  trial  and  error  that  required  only 
two  iterations  for  large  drops  and  no  more  than  six  for  small  drops  in  order  to  converge  to  a  solution 
of  the  boundary  value  problem. 


The  constants  Kj  and  K2  group  the  aerodynamic  and  gravitational  parameters  respectively. 
The  data  from  References  3  and  5  for  values  of  C  are  given  in  Table  1  for  various  values  of  Re . 

Equations  (8)  to  (11)  are  nonlinear  ordinary  differential  equations  for  the  drop  co-ordinates 
x  and  z.  The  local  flow  velocities  u„  and  w,  are  not  functions  of  time  in  steady  flight.  They  describe 
the  spatial  variation  of  file  local  flow  effects  from  the  wing  and  PMS  probe  canister.  Solutions  to  the 
above  system  must  be  performed  numerically  because  of  the  variation  of  C  with  Re  and  the  variation 
of  u,  and  wt  with  x  and  z. 

The  initial  values  of  the  droplet  velocity  components  are 


r  v. 


The  initial  position  of  the  droplet  must  be  chosen  such  that  its  final  position  is  within  the 
probe  sample  volume.  Thus,  the  boundary  value  problem  is  posed: 

find  the  initial  position  of  the  drop  with  a  given  initial  velocity  such  that 
the  final  position  condition  is  satisfied  at  an  unspecified  velocity. 


2.1  Initial  Portion 

Aa  stated  before,  the  initial  position  of  the  drop  must  be  found  so  that  its  final  position  is 
within  a  specified  distance  of  the  top/bottom  of  the  sample  volume.  The  initial  x  co-ordinate  of  the 
drop  was  taken  to  be  31  quarter-chords  upstream.  The  initial  z  co-ordinate  was  varied  to  solve  the 
boundary  value  problem.  The  initial  x  was  chosen  to  be  far  enough  away  from  the  wing  so  that  the 
flow  effects  due  to  the  wing  were  negligible.  To  prove  this,  one  set  of  calculations  with  an  initial  x  of 
-21  was  computed  and  the  results  showed  no  discernible  difference  from  those  obtained  using  an 
initial  x  of -31.  Similar  conclusions  were  drawn  when  an  initial  x  of -41  was  used.  All  the  results 
presented  here  used  an  initial  x  of -31. 


2.2  Find  Position 

A  successful  solution  to  the  boundary  value  problem  occurs  when  the  final  position  of  the 
droplet  is  “sufficiently”  dose  to  the  desired  value.  In  this  work,  a  computational  box  in  the  x,  z  plane 
was  placed  around  the  final  point.  The  dimensions  of  the  box  relative  to  the  central,  final  point  were 
±  0.01  quarter-chords  in  the  x  direction,  and  ±  35 /a  or  ±  2  X  Iff*4 ,  which  ever  was  bigger,  in  the  z 
direction. 


All  of  the  computations  were  performed  for  the  case  of  a  PMS  2D-C  probe  mounted  on  the 
short  pylon  at  the  inboard  location  along  the  span.  The  dimensions  and  sample  volume  location  are 
given  in  Table  3.  The  probe  counts  droplets  that  are  within  its  depth  of  field.  Hence,  both  the  top  and 
bottom  of  tiie  sample  volume  had  to  be  considered  as  final  positions  for  the  drops. 


3.0  THE  LOCAL  FLOW  VELOCITIES 

The  local  flow  velocities  u,  and  w,  are  functions  of  position.  These  velocities  were  calculated 
by  modeling  the  flow  about  the  Twin  Otter  wing  using  Joukowski  aerofoil  therory  (Ref.  2)  and  in¬ 
cluding  canister  effects.  The  validity  of  the  model  was  established  by  measuring  the  flow  velocity  at 
the  four  locations  where  probes  are  carried  under  the  Twin  Otter  wings. 

A  brief  summary  of  the  results  of  Reference  2  is  included.  The  Twin  Otter  wing  was  taken 
to  have  a  rectangular  planform,  a  symmetrical  section  of  16%  thickness,  an  elliptical  span  wise  load 
distribution  and  no  twist.  The  simplest  Joukowski  transformation  was  employed  which  caused  the 
model  of  the  real  wing  to  have  a  cusped  trailing  edge  instead  of  having  the  correct  trailing  edge  angle. 
The  forward  half  of  the  wing  section  was  well  represented  as  a  Joukowski  aerofoil. 

The  main  conclusion  of  that  report  was  that  the  Joukowski  aerofoil  theory  (including 
canister  effects)  could  be  used  to  obtain  u,  and  wa  anywhere  ahead  of  the  wing  quarter-chord  point 
with  an  accuracy  bettor  than  97%,  for  the  probe  location  of  Table  3,  in  the  range  of  aircraft  lift  co¬ 
efficients  (CL  )  applicable  to  cloud  penetrations.  The  lift  coefficient  is: 


CL  -  L/(l/2  pV2S) 

For  straight  and  level  flight,  the  lift  L  equals  the  aircraft  weight. 


A  point  in  the  wing  plane  is  located  by  the  complex  co-ordinate  q 
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and  the  corresponding  complex  co-ordinate  f  in  the  cylinder  plane  is  located  by  the  conformal  trans¬ 
formation 


f  -  (tj/2)  ±  V(r?/2)2  -  1 


(22) 


The  sign  is  chosen  so  that  f  is  outside  the  cylinder.  The  velocity  in  the  wing  plane  is 
calculated  most  easily  by  the  relation: 


u,  "  1W, 


[ 


exp(-  ia) 


(1  +  c/a)2 
tf  +  e/a)2 


exp(ia)  +  2i 


(1  +  e/a) 
(£  +  e/a) 


(23) 


where  the  overbar  on  the  x  component  of  velocity  signifies  that  only  flow  due  to  the  wing  is  described 
by  Equation  (23). 

The  thickness  parameter  is  e  and  the  wing  angle  of  attack  is  a.  A  method  of  obtaining  e  for 
a  given  wing  thickness  is  described  in  Reference  2. 

The  effects  of  the  canister  were  included  in  a  simple,  approximate  fashion.  The  discussion 
in  Reference  2  on  the  modeling  of  the  flow  ahead  of  the  hemispherical  end  of  the  canister  can  be 
summarized.  The  flow  due  to  the  wing  in  the  x  direction,  u,,  is  modified  by  the  flow  about  a  sphere. 
The  following  relation  is  used  for  the  final  estimate  of  the  local  velocity  u,,  because  the  canister  is 
considered  to  be  “small’*  and  imbedded  in  the  flow  about  the  wing: 

“.  -*.(i-<*A>5) 

where  r(  is  the  radius  of  the  spherical  end  of  the  canister  and  x,  is  the  distance  of  the  drop  from  the 
carter  of  the  sphere.  The  small  effects  of  the  sphere  on  wa  are  neglected  and  u,  is  taken  to  be  along 
the  axis  of  the  canister.  Both  r,  and  x,  are  in  terms  of  the  fundamental  length  ‘a’. 


4.0  CLOUD  PHYSICS  IMPLICATIONS 

There  are  two  main  implications  of  aircraft-induced  flow  effects  on  cloud  physics  results. 
The  first  is  that  the  images  of  drops  formed  by  PMS  probes  are  distorted  unless  the  drops  pass  through 
the  probe  sample  volume  at  the  true  airspeed  of  the  aircraft.  The  scanning  rate  of  the  laser  in  the  PMS 
probe  is  adjusted  to  give  the  correct  images  to  drops  under  the  assumption  that  they  are  moving  at  the 
aircraft  speed  and  not  at  a  different  speed  because  of  local  flow  effects.  For  example,  if  the  drop 
velocity  is  80%  of  aircraft  true  airspeed  as  it  passes  through  an  imaging  probe,  then  its  dimension  in 
the  x  direction  will  be  distorted  by  the  ratio  of  five  to  four  relative  to  the  spanwise  dimension.  The 
drop  will  appear  elongated  in  the  x  direction. 

Since  the  velocity  of  a  drop  that  passed  through  the  bottom  of  the  sample  area  is  in  general 
different  from  the  velocity  of  one  of  the  same  size  that  passed  through  the  top,  then  some  rule  must 
be  developed  for  chosing  the  corrective  factor  to  be  applied  to  images.  Both  drops  will  be  in  focus  and 
will  be  imaged  but  there  is  no  way  to  tell  which  image  goes  with  a  “top  drop”  or  a  “bottom  drop”.  It 
is  suggested  that  the  average  of  the  velocities  be  used.  Acceptable  accuracy  is  possible  only  if  the  “top” 
and  “bottom"  velocities  are  not  too  different. 
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The  second  effect  is  that  measurements  of  drop  concentrations  (number  per  unit  volume) 
are  distorted  because  drops  of  different  sizes  as  counted  by  the  probes  originate  from  cloud  volumes 
of  different  sizes  and  different  locations.  It  is  usually  assumed  that  drops  move  on  straight  lines  so 
that  the  volume  swept  out  by  a  probe  is  determined  only  by  the  frontal  dimensions  of  its  sample  area 
and  the  aircraft  speed.  However,  small  drops  tend  to  follow  the  local  flow  streamlines  and  if  the 
streamlines  at  infinity  expand  (or  contract)  relative  to  their  spacing  at  the  probe,  then  the  volume  used 
in  the  concentration  calculation  is  incorrect.  Large  drops  ignore  local  flows  and  pass  through  the  probe 
sample  area  as  though  they  had  moved  on  straight  lines  (almost)  that  remain  parallel.  In  this  case,  the 
assumed  volume  used  in  the  concentration  calculations  is  the  correct  one.  Concentration  measurements 
become  a  function  of  drop  size  because  trajectories  of  different  sized  drops  have  different  starting 
points. 


5.0  RESULTS  FOR  VELOCITY 

Drop  trajectories  that  solved  the  boundary  value  problem  wore  computed  for  one  probe 
location.  Trajectories  passing  through  the  top  and  bottom  of  the  sample  area  were  studied  and  the 
drop  velocity  as  it  passed  through  those  positions  was  computed.  At  the  same  time,  the  ratio  of  the 
difference  in  the  initial  z  co-ordinate  to  the  difference  in  the  final  z  co-ordinate  of  trajectories  through 
the  top  and  bottom  of  the  sample  areas  was  calculated.  This  calculation  is  pertinent  to  analysis  of 
concentration  errors  while  the  former  is  related  to  image  distortions. 

The  velocity  ratio  is 


VR  ”  yj  u2  +  w2 


If  VR  is  less  than  one,  then  the  droplet  is  moving  dower  than  the  aircraft  speed. 

Figures  4, 5  and  6  show  three  typical  results  for  the  droplet  velocity  ratio  at  the  top  and 
bottom  of  the  sample  area  as  a  function  of  drop  size.  Each  figure  is  applicable  to  one  set  of  flight 
conditions  (aircraft  lift  coefficient). 

In  all  cases,  the  velocity  ratio  of  the  largest  drop  tends  to  one  while  the  ratio  for  the  smallest 
drop  approaches  that  of  the  local  airflow.  Both  of  these  results  were  expected.  As  aircraft  lift  coeffi¬ 
cient  increases,  the  velocity  ratio  for  the  small  drops  decreases  and  the  local  minimum  for  moderately 
sized  drops  becomes  less  pronounced.  The  minimum  occurs  at  a  smaller  drop  size  for  increasing  CL 
but  it  usually  occurs  for  drop  diameters  between  40  and  80p.  The  exact  nature  of  the  curve  in  this 
region  depends  on  the  integrated  effects  of  the  velocity  field  experienced  by  the  drop  as  it  approaches 
the  probe. 


The  cause  of  the  minimum  in  Figures  4  to  6  is  the  inertia  of  the  drop.  As  the  drop  ap¬ 
proaches  the  PMS  probe,  there  is  a  decreasing  magnitude  of  the  velocity  field  due  to  the  wing.  The 
inertia  of  the  droplet  does  not  let  it  slow  as  quickly  as  the  local  flow.  Then,  the  velocity  field  starts  to 
accelerate  and  the  drop  lags  this  trend  because  of  its  inertia.  Hence,  it  goes  through  the  sample  volume 
at  a  velocity  less  than  that  of  the  local  flow.  The  above  discussion  is  not  true  for  drops  so  small  that 
they  essentially  follow  local  flow,  nor  for  drops  so  large  that  the  local  flow  is  of  no  consequence. 


In  all,  nine  values  of  CL  were  chosen  for  computation  of  the  velocity  ratio  at  the  top  and 
bottom  points  of  the  probe  sample  volumes.  The  range  of  variables  was  V  ■  60,  70,  80  m/s  and 
h  ■  6,000, 10,000  and  14,000  ft.  and  the  aircraft  weight  was  chosen  as  10,000  lbs  in  all  cases.  The 
nine  values  of  the  corresponding  CL  are  shown  in  Table  4.  Results  for  only  three  of  these  were  shown 
in  Figures  4  to  6  but  th«y  are  typical  of  the  other  six. 
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For  a  constant  diop  size,  the  avenge  velocity  ratio  can  be  plotted  versus  CL  for  drops  that 
PM*  through  the  sample  volume.  This  procedure  was  adopted  for  a  few  drop  sizes  and  the  results  are 
shown  in  Figures  7  and  8.  Figure  8  applies  to  small  drops  and  shows  the  minimum  discussed  above. 

A  least  sqm  stm  fit  to  the  nine  values  was  made  and  the  straight  line  is  the  result  of  the  procedure. 

Close  examination  of  file  figures  reveals  that  the  deviation  of  each  point  from  the  line  is  small,  but 
that  H«t«  at  each  speed  really  fall  into  separate  populations.  The  aircraft  CL  does  not  supply  exact 
correlation  to  the  data.  The  reason  is  that  the  local  flow  velocity  is  a  linear  function  of  CL  but  the 
droplet  drag  depends  on  atmospheric  viscosity  as  well.  However,  the  correlation  with  lift  coefficient 
is  considered  to  be  adequate. 

Figure  9  is  a  plot  of  all  the  velocity  ratios  versus  CL  and  drop  size  for  the  average  of  the  top 
and  bottom  velocities  at  the  probe  sample  volume.  This  plot  masks  the  slight  variation  of  drop  velocity 
with  speed  and  altitude.  The  dotted  line  in  the  CL  -diameter  plane  shows  where  the  straight  lines  for 
smaller  drops  would  pierce  the  plane  at  high  CL.  This  figure  shows  how  the  local  minimum  moves  to 
smaller  drop  sizes  as  CL  increases.  The  results  are  extrapolated  to  a  large  value  of  CL  (1.2).  The  local 
flow  velocity  model  is  not  very  reliable  in  this  area  because  the  wing  is  approaching  a  stalled  condition. 

Figure  10  shows  the  intercepts  and  slopes  of  the  straight  lines  drawn  in  Figure  9.  The 
smoothness  of  the  curves  show  that  interpolation  between  the  data  points  is  valid  and  that  limited 
extrapolation  to  drops  larger  than  1  mm  in  diameter  is  permissible. 


6.0  RESULTS  FOR  TRAJECTORY  EXPANSION 

The  trajectory  expansion  ratio  ER  is  defined  as: 

ER  *  Azj/Azf 


where  Az  is  the  positive  difference  in  position  between  “top”  and  “bottom"  trajectories  previously 
described  and  the  subscripts  “i”  and  ‘Y’  refer  to  initial  and  final  position  respectively.  If  ER  is  less 
than  one,  there  has  been  “expansion”  from  the  initial  to  the  final  position. 

Figures  11, 12,  and  13  show  the  expansion  ratio  for  the  same  values  of  CL  that  applied  to 
the  previously  described  velocity  ratio  results  (Figs.  4, 5  and  6).  There  is  more  scatter  in  the  data  be¬ 
cause  both  the  numerator  and  denominator  in  the  expression  for  ER  are  obtained  from  the  difference 
of  two  small  numbers  of  the  same  order  of  magnitude.  Also,  a  solution  is  accepted  if  a  drop  is  within 
the  computational  box  which  admits  a  tolerance  on  Xf  and  if.  However,  the  smoothness  of  the  data  is 
sufficient  to  allow  the  ER  as  shown  in  file  figures  to  be  used  to  correct  cloud  measurements.  These 
results  are  restricted  to  the  case  of  spherical  drops. 

At  low  CL ,  Figure  11,  there  is  almost  no  drop  concentration  distortion  due  to  local  flow 
effects  because  ER  is  very  nearly  equal  to  one  for  all  drop  sizes.  As  CL  increases,  (Figs.  12  and  13)  ER 
decreases  for  small  drops  and  remains  near  one  for  large  drops.  Underestimation  of  drop  concentra¬ 
tions  is  the  general  result. 

The  data  from  the  above  three  figures  can  be  plotted  versus  CL  for  a  given  drop  size 
(Figs.  14  and  16).  Figure  16  applies  to  smaller  drops.  A  least-squares  straight  line  fit  to  the  nine  data 
points  is  shown.  Very  small  deviations  of  data  from  the  line  are  evident.  The  results  are  qualitatively 
the  same  as  those  for  the  velocity  ratio  but  there  are  small  quantitative  differences.  If  desired,  a  plot 
similar  to  Figure  9  could  be  drawn  to  suit  the  expansion  ratio. 

Figure  16  shows  the  slopes  and  intercepts  of  the  straight  lines  obtained  in  Figures  14  and  16 
when  plotted  versus  drop  diameter.  Again  there  is  more  scatter  when  compared  to  the  similar  plot  of 
velocity  ratio  parameters  (Fig.  10).  The  trends  and  shapes  of  the  curves  on  Figure  16  are  about  the 
aame  as  those  on  Figure  10. 


7.0  CORRECTIONS  TO  CLOUD  PHYSICS  MEASUREMENTS 


Drops  passing  through  the  sample  volume  at  speeds  other  than  the  aircraft  speed  have 
images  that  are  distorted.  In  general,  the  images  of  all  drops  will  need  correction.  The  drop  dimension 
in  the  x  direction  must  be  multiplied  by  the  velocity  ratio  appropriate  to  the  drop  diameter  and  air¬ 
craft  lift  coefficient.  The  diameter  5  is  obtained  from  the  drop  dimension  perpendicular  to  the  x,  z 
plane,  which  is  given  by  the  photo-diode  array.  Sufficient  data  to  compute  aircraft  CL  must  be  re¬ 
corded  at  an  adequately  high  sampling  rate  so  that  the  x  dimension  of  each  drop  can  be  scaled  in  the 
software  at  the  CL  appropriate  to  the  time  of  diameter  measurement. 


Let: 


VR(5,Cl)  =*  A(6)  +  B(S)Cl 


where  A(5)  and  B(5)  are  the  slope  and  intercept  of  a  straight  line  as  illustrated  in  Figures  7  and  8. 
These  two  parameters  are  plotted  in  Figure  10.  Interpolation  functions  could  be  determined  empiri¬ 
cally  for  A(5)  and  B(5)  so  that  if  one  knew  6  and  CL ,  the  correct  velocity  ratio  could  be  calculated. 
Conversely,  the  data  of  Figure  10  could  be  used  to  form  a  table  of  values  for  various  values  of  8.  This 
table  could  be  stored  in  the  computer  and  used  for  on-board  analysis  or  for  later  analysis  on  the 
ground.  At  each  value  of  8,  interpolation  between  entries  in  the  table  would  lead  to  an  approximation 
to  the  correct  values  of  A  and  B.  The  coarseness  of  the  approximation  would  depend  on  the  size  of 
the  interval  between  table  entries.  Perhaps  the  time  available  in  the  computation  cycle  would  determine 
these  intervals  and  this  could  dictate  the  accuracy  of  image  correction  caused  by  drop  velocity  effects. 

Measurements  of  drop  concentrations  can  be  corrected  in  the  same  way  by  using  A(6)  and 
B(8)  from  Figure  16.  Division  of  concentration  data  by  the  expansion  ratio  ER  allows  the  effects  of 
local  flow  to  be  accounted  for.  Without  this  correction,  concentration  data  are  under-estimated  if  ER 
is  less  than  one  and  over-estimated  if  ER  is  greater  than  one.  Small  drops  at  low  CL  exhibit  the  latter 
error  while  most  drops  experience  the  former  error. 

The  numerical  results  of  the  calculations  presented  in  this  report  apply  to  one  PMS  probe 
installation  on  the  Twin  Otter:  the  inboard  spanwise  position  with  the  short  pylon.  These  calculations 
can  be  repeated  for  the  other  three  possible  PMS  mounting  locations  in  order  to  obtain  corrections  to 
measurements  at  those  locations.  The  purposes  of  this  report  have  been  fulfilled  without  doing  these 
calculations.  In  any  event,  the  computer  programs  and  data  reduction  logic  have  been  developed.  The 
calculations  could  be  routinely  performed  and  would  not  show  large  differences  from  the  results 
quoted  here  because  the  flow  is  not  greatly  different  at  those  points. 

In  the  range  of  CL  between  0.3  and  0.76,  the  concentration  measurements  of  small  drops 
are  at  most  18%  under-estimated  at  present.  The  concentration  measurements  of  large  drops  are 
approximately  correct.  The  same  is  true  for  images,  which  are  at  worst  distorted  by  only  20%  in  the 
axis  ratio  determination.  These  errors  are  not  large  but  they  can  be  nearly  removed  by  applying  the 
methods  described  here  using  Figures  10  and  16  to  obtain  the  slope  and  intercept  for  each  drop  size. 
Application  of  these  parameters  to  the  measured  aircraft  lift  coefficient  should  lead  to  image  and  con¬ 
centration  measurements  correct  to  a  few  percent.  The  theoretical  local  flow  velocity  is  more  nearly 
equal  to  the  measured  value  at  two  other  probe  locations  (Ref.  2)  so  one  would  expect  even  better 
corrections  for  probes  mounted  here. 

The  use  of  CL  based  on  aircraft  weight  to  correct  cloud  physics  measurements  can  be 
criticized  in  so  far  as  die  local  flow  is  not  completely  steady  because  of  aircraft  response  to  gusty  air 
within  a  cloud.  However,  some  of  the  in-flight  measurements  described  in  Reference  2  were  taken 
under  “gusty”  conditions  and  the  effect  of  this  unsteady  flow  on  the  correlations  presented  there 
could  not  be  noticed.  It  is  stated  without  proof  that  unless  the  aircraft  motion  becomes  very  unsteady, 
the  methods  presented  here  will  still  apply. 
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Errors  in  imaging  and  concentration  measurements  can  arise  from  many  sources  and  only 
those  resulting  from  aircraft  flow  effects  have  been  discussed.  It  is  beyond  the  scope  of  this  report  to 
attempt  a  general  treatment  of  all  types  of  errors. 


8.0  CONCLUSIONS 

Cloud  physics  measurements  by  the  Twin  Otter  aircraft  are  affected  to  some  degree  by  local 
flow  effects  from  the  wing  and  from  the  probe  canister.  Droplet  images  and  concentrations  are  dis¬ 
torted  in  the  worst  case  by  no  more  than  20%  for  an  aircraft  lift  coefficient  of  0.75  and  by  only  a  few 
percent  at  lower  values  of  CL .  Typically,  doud  penetrations  take  place  at  a  CL  of  between  0.4  and 
0.6.  Corrections  for  a  Partide  Measuring  Systems  2D-C  probe  mounted  at  the  inboard  location  on  the 
short  pylon  have  been  given  which  are  based  on  an  experimentally  verified  model  for  the  local  flow 
vdodty. 
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T ABLE  2  I 
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ATMOSPHERIC  PROPERTIES 
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v  ft2  /sec 
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cft/see 

pduas/lt3 

P  mb 

Sea  level  (SL) 

1.576  X  10"4 

518.69 

1116.4 

1013.2 

6000 

1.869  X  10“4 

1 

14000 

2.234  X  10'4 

TABLES 

DIMENSIONS  AND  LOCATION  OF  SAMPLE  VOLUME 


Sample  volume  co-ordinates  (in  quarter  chords) 


Position 

z 

a 

Center 

-1.595 

-0.58 

Top 

-1.595 

-0.518 

Bottom 

-1.595 

-0.642 

Canister  nose  center 

-1.130 

-0.58 

Wing  Chord:  78  inches 

Probe:  PMS  2D-C  mounted  at  61.8%  of  semi-sDan 
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FIG.  4:  VELOCITY  AT  PROBE 


FIG.  7:  VELOCITY  AT  PROBE 
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FIG.  8:  VELOCITY  AT  PROBE 
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FIG.  10:  VELOCITY  RATIO  PARAMETERS 
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FIG.  11:  TRAJECTORY  EXPANSION 
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FIG.  IB:  TRAJECTORY  EXPANSION 
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